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parameter specifying the amount of vortex rollup 
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local longitudinal loading obtained by spanwlse integration of 
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number of wake corner points 

free-stream Mach number 
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cot A 

radius of the circular arc in the transformed plane 

radius of the circle in the transformed plane containing the vor- 
ticity in the spiral region in the Mangler-Smlth theory 
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resultai.t velocity field as computed by the Potential Flow 
AnaT'Sls program (POTFAN) 
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rectangular coordinate system fixed in the canard-wing plane, 

X increases in the stream direction, y horizontally to 
starbot.rd, and z upward 

position vector relative to the x, y, z-coordinate system 
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angle of attack, deg 
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t>aran»cers characterizing the vortlclty distribution 
singularity strength 

C (upper surface) - C (lower surface) 
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element of length along panel boundary 
pressure difference 

half angle subtended by the finite circular arc representing the 
free sheet In the transformed plane 

angular coordinate on this circular arc 
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measure of the effective vortlclty within the spiral region of 
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VORTEX EFFECTS FOR CANARD-WING CONFIGURATIONS AT HIGH ANGLES 


OF ATTACK IN SUBSONIC FLOW 
B. M. E. de Silva* and R. T. Medan 
Amea Research Center 


SIFMMARY 


A fully three-dimenHlonal subsonic panel nM*thod that ran handle arbi- 
trary shed vortex wakes is used to compute the nonlinear forces and moments 
on simple canard-wing configurations. The lifting surfaces and wakes are 
represented by doublet panels. The Mangler-Smlth theory is used to provide 
an initial estimate for the vortex sheet shed from the leading edge. The 
tralling-edge wake and the leading-edge wake downstream of the trailing edge 
are assumed to be straight and leave the trailing edge at an angle of u/2. 

Results indicate good agreement with experimental data up to 40* angle 
of attack. 


INTRODUCTION 


The object of the research described herein is to contribute to the 
development of a method for predicting the subsonic nonlinear forces on com- 
plete canard-wing aircraft configurations with leading-edge separation and 
vortex rollup. In particular, a simplified technique for estimating the 
shape and vortlcity distribution of vortex sheets generated by simi)le canard- 
wing configurations is described. Although the simplified technique described 
herein will be used to generate an initial estimate for a more precise itera- 
tive technique that will be developed at a later stage, it is sufficiently 
accurate In certain cases to justify its use without subsequent Improvement 
of the wake shapes. 


METHODOLOGY 


Leading-Edge Vortex Sheet 

The method used for estimating the shape of the wake is based on the 
Mangler-Smlth theory (refs. 1-3). It assumes conical, slender body theory in 
the crossflow plane giving a two-dimensional Laplacian as the governing equa- 
tion. The boundary conditions are the vanishing of the normal component of 
velocity on the wing and vortex sheet and continuity of pressure across the 
sheet. 
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The analysis Is fornulated in a transforned plane on which the wake 
sheet emanating from the leading edge called the free sheet Is approximated 
by a circular arc of radius r^ and arc length 2 rj6 (fig. 1(a)). This arc 
then connects continuously at point C* to a small circle of radius 
representing the so-called fed sheet and containing the vorticlty p. The 
vorticity distribution Y on the free sheet is assumed to have the forr 

Y(e') - Yo + Yi cos 0' - Yj sin 0’ (0 < 0' < 20) (1) 

Thus the shape of the sheet is described by r^, r 2 t ^ while its vorticlty 
is represented by u'» Yq» Yi, Y 2 . These seven parameters are determined by 
enforcing the boundary conditions at isolated locations, together with the 
requirements of self-consistency and compatibility. This results in a system 
of nonlinear equations which are solved through a least squares minimization 
(ref. 4) with simple inequality constraints of the form 

0 < t2 < rj ; p > 0 ; 0 < 0 < I (2) 


Initial Wake Estimate 

The Mangler-Smith theory gives an initial estimate of the shape of the 
separated sheet with reasonable accuracy. To compute the flow field, this 
sheet is discretized so that the separated wake is essentially modeled by a 
collection of concentrated line vortices, which are described using vortex 
quadrilaterals. This vortex description is then combined with a panel code 
in which the wakes are represented by doublet panels which are either constant 
in strength or have a quadratic variation. (The equivalence between doublet 
panel edges and line vortices alluded to here is discussed in the subsequent 
section on panel methods.) 

To align the vortex trajectories with the streamlines, the corner points 
along the leading edge of each wake panel are obtained by integrating the 
velocity field. Specifically, the flow tangency condition along the leading- 
edge wake is given by 


dX 

^ - V - (V • N) N (3) 

where N is the unit normal to the vortex sheet and ^ is the velocity field 
as computed by the Mangler-Smith theory. This consists of the free-stream 
velocity together with the velocities induced by the vortex sheet and its 
image. The corner points ; j ■ 1.2, . . .,J) are obtained by a 

numerical integration of equation (3) based on a modified fourth-order 
Runge-Kutta technique. Each shed wake vortex is represented by J-1 seg- 
ments where the value of J is such that the wake vortices can reach the 
trailing-edge plane. The traillng-edge wake and the leading-edge wake 
downstream of the trailing edge are assumed straight xeaving the trailing 
edge at an angle of a/2. 
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Panel Methods 


It ir emphasized that the technique described in the preceding sections 
is used only for estinating the leading-edge wake shape and not for determin- 
ing its vorticity distribution or effect on pressures. Given the wake shape, 
the latter is accomplished by a fully three-dimensional technique that may be 
described as a panel singularity method based on the linear differential 
equation 


(1 - M*) 0 (4) 

XX yy zz 

For purposes of comparison, two panel codes were used. The first method 
uses constant-strength doublet panels and is referred to as the Potential Flow 
Analysis Program or POTFAN (refs. 5-10). The second method uses quadrat Ically 
varying doublet panels and is referred to as the Advanced Panel Code or APC 
(refs. 11 and 12). The APC is still in the process of develo(xnent , conse- 
quently its range of applicability was limited. Therefore, the primary 
emphasis in this paper is centered on results obtained from POTFAN. 

In the first method, POTFAN, the lifting surfaces and wakes aie repre- 

sented by constant doublet panels. Away from the lifting surface edges, the 
panels are quadrilaterals with all four edges bound in the surface. Panels 

on the edges have two or three straight edges bound to the surface with two 

other edges in the flow field. The edges in the flow field are comprised of 
many short, straight line vortex segments plus a final straight segment that 
goes off to infinity. Figure 1(b) is a schematic of the doublet panels used 
to represent a lifting surface with flow separation ou all sides. 

The kinematic and dynamical boundary conditions on the wake are satisfied 
if the panel edges in the flow field are streamlines. It is currently assumed, 
that the flow-field panel edges are sufficiently close to streamlines. Given 
this assumption, the strengths of each of the panels may be determined by 
enforcing the v * ii ■ 0 boundary condition at one point inside each panel. 

If there are M panels on the lifting surface, this boundary condition 
results in a system of M linear equations in M unknowns. Once these 
equations are solved for the doublet strength for each panel, the pressure 
distribution can be calculated using equation (5). 

AE - pTjf * At (5) 

The forces and moments are obtained by integrating the resulting pressure 
distribution. 

In calculating the velocity induced by any of the constant doublet 
panels, it is recognized that such a panel is equivalent to a concentrated 
line vortex situated at the edge of the panel and therefore, the Biot-Savavt 
law is available for calculating the panel-induced velocities. Thus the wakes 
may be thought of as being represented by concentrated line vortices. If the 
center of these line vortices pass close to lifting surfaces, the unmodified 
Blot-Savavt law may yield unrealistically high induced velocities. Therefore, 
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the concept of a vortex core ia used. Within each core the velocity within 
the core varies linearly with the core radius (l.e., the cores are Ranklne 
cores). The cores were taken to be constant In diameter and the diameters 
are generally chosen to be the sane as the panel widths at the leading edge. 

The second method Is a higher order panel code (refs. 11-12), which uses 
linearly varying source and or quadrat Ically varying doublet distributions 
over each of the panels. This method Is applied In essentially the same 
manner as the first method. However, with the quadrat Ically varying doublet 
distributions on the wake panels, the wake Is a continuous vortex sheet and 
there Is no need to Introduce the concept of vortex cores Into the method. 


Configuration Modeling 

To evaluate this technique It was applied to both Isolated wings and a 
canard-wing combination. One Isolated wing was Investigated to evaluate the 
details of the methods which Included the load distributions, the flow-field 
velocities above the wing, and the overall forces and moments. The aspect 
ratio and leading-edge sweep of the two wings were nominally 1.0 (1.15 and 
1.07) and 75* (74* and 75*), respectively. The canard-wing combination was 
Investigated to assess the ability of the method to predlc' the vortex* 

Induced forces on these types of configurations. The study considered the 
effects of the canard wake for a simple canard-wing configuration of trape- 
zoidal planform with the canard and wing In the same plane. It Is assumed 
that inboard of the canard tip the t railing-canard wake produces a downwash 
over the wing which counter.icts the angle of attack, thus decreasing the 
leading-edge suction and eliminating leading -edge separation frmn this portion 
of the wing. Outboard of the canard tip there Is an upwash field which 
Increases the wing local angle of attack, thus flow separation Is assumed 
outboard of the canard tip. The model used la shown In figure 2. The experi- 
mental results are contained In reference 13. 


RESULTS 


The results are displayed in two parts. Figures 3-9 show the results 
for flat delta wings with sharp leading edges, while figures 10-15 show the 
results for the canard-wing combination of figure 2. Figure 3 shows a spiral 
vortex for a deita wing of ^ ■ 0.7 at an angle of attack of 40*. The 
assumption of conical flow In the Mangler-Smith theory implies at each 
x-statlon, the leading-edge vortex trajectories have the property that the 
family of curves z/x vs y/x are coincident. Figure 4 shows different views 
of the vortex trajectories. As previously noted, the leading-edge wake 
vortices are assumed to be straight downstream of the trailing edge leaving 
the wing at an ..ngle of a/2 and extending to downstream Infinity. In 
addition, the shed wake from the apex goes along the vortex core. 

Figures 5-8 show the normal force characteristics for a wing with an 
Al ■•1.15. The nominal flow model consisted of 256 wing panels for POTFAN 
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and 225 for Che APC. Figure 5 comparea the lift coefficient predicted by 
POTFAN and APC with the Polhaamua analogy (ref. 14) and experiment ( ref. 15) 
up Co 40* angle of attack. The agreement among the panel codea and experi- 
mental reaulta ia very good. It la alao aeen that Che Polhamoua auction 
analogy ia quite good for predicting the lift. 


Figure 6 comparea the longitudinal load diatribution for a > 10*. 20*, 
and 30* with experiment. The local normal force ia obtained by integrating 
the apanwiae preaaurea at each chordwiae ataclon. 



( 6 ) 


Thia integration la carried out numerically by sumring the preaaure differ- 
encea over the panela in the apanwiae direction. Again, the predicted values 
agree well with experiment. 


Figure 7 explores the effect of increasing the number of panels indi- 
cating that Che method is convergent insomuch as Che agreement with experi- 
ment improves as the number of prnels increases. 


Since Che trail ing-edge Kutta condition has not been strictly enforced 
(ref. 16), the center of pressure is not as far forward as the experimental 
position as shown in figure 8. In conical flow the center of pressure is at 
the centroid of Che wing and figure 8 also shows the three-dimensional panel 
singularity field to differ substantially from conical flow. 

The APC has a current upper limit of 1000 panels including wake panels. 
The total number of panels used were 1496 for POTFAN and 990 for the APC. 

To get better agreement with pressure and moment data, the APC would require 
more panels Chan currently allowed. This limitation also precludes the APC 
operating at a < 20*, as Che Mangler-SmiCh technique requires an increased 
number of wake panels to adequately represent the rollup as Che angle of 
attack decreases. POTFAN, on Che other hand, can operate down Co about 5-10*. 

Finally, figures 9(a), (b) , and (c) compare the spanwise variation of the 
velocity components predicted by POTFAN with measurements from reference 17 
at flow-field points above e. wing of A ■ 1.07 at a ■ 29*. At the time of 
this study, the pilot version of the APC did not allow for velocities Co be 
computed at off-body points so no comparable results are available from Che 
APC. The predicted core location in the crossflow plane is y/s •• 0.693 and 
z/s ■ 0.416, while the experimental value is y/s ■ 0.710 and z/s ■ 0.413. 

The axial velocity, figure 9(a), exhibits steep gradients as the core is 
approached from either side and a leveling off near the center line y •• 0. 
Near Che surface (z/s ■ 0.102), the velocity exhibits some small fluctuations 
due probably to boundary-layer effects. Near the surface the spanwise 
scrubbing action of the vortex is shown in figure 9(b), with the peak lying 
directly under the vortex core. AC the level of the vortex core, there is a 
perceptible discrepancy between predicted and experimental values. More 
accurate calculations (ref. 18) predict Che cross section of the fed sheet to 
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be flattened with the long axle tilted top Inward. Thla dlatortlon from the 
roughly circular aplral predicted by the preaent formulation would at leaat 
account for the change In direction of the spanwlae component as the core la 
pasacd. The vertical velocity, figure 9(c), again shows the effects of strong 
vortex flow. At the core height, z/s ■ 0.413, downflow occurs Inboard of the 
core and upfiow outboard. Near the surface, the downflow over the inboard 
section goes to zero at the span-wise location of the core. The flow then 
gets directed upward /out board of the core. 

The programs were run on a CDC 7600 with typical run times of about 
3.6 min for POTFAN for a wing with around 1500 panels (wing and wake panels). 
Corresponding APC runs with about 990 panels required 6.8 min. 

Figures 10-15 compare the predicted and measured results for the canard- 
wing configuration shotm In figure 2. Figure 10 shows the leading-edge wake 
for the canard and wing. Note the wing leading-edge wake beginning outboard 
of the canard. It should also be noted that the Mangler-Smlth theory gives 
an adequate description of the rollup process only If the parameter 

a - 0.0175 a tan A (7) 

Is sufficiently large. This Implies that unless the leading edge Is highly 
swept, a must be large. For this reason, the rollup can be effectively 
modeled only for a > 24*. 

Figure 11 shows the canard and total lift coefficients up to a • 40*. 
Below a ■ 24*, the rollup from the leading edges Is omitted. For purposes 
of simplicity. In this preliminary Investigation, the body Is not Included 
In the rollup calculations. Figure 12 shows the drag polars. Good agreement 
with test data In figures 11 and 12 Indicates that the rollup effects are 
small for a < 24*. There is a significant discrepancy between the predicted 
and measured pitching-moment characteristics as illustrated in figure 13. 

The pitching moments are probably more sen Itlve to the absence of a fuselage, 
thus accounting for some of this discrepancy. 

Figures 14 and 15 give the pressure distributions Inboard and outboard of 
the canard for o ■ 28*. The vortex cores for the can.ird are situated near 
its leading edge so that the pressure coefficients become quite large near 
the canard leading edge and drop rapidly as the midchord Is approached 
(fig. 14). The pressures on the wing are substantially lower In the vicinity 
of Its leading edge, gradually leveling off and tending toward zero as the 
trailing edge Is approached. Outboard of the canard (fig. 15), the pressure 
does not tend to zero in the neighborhood of the trailing edge, due to the 
Kutta condition not being satisfied. 

A typical POTFAN run for the wing-canard case needed a total of about 
5700 panels with conventional paneling of the type shown in figure 2. Of 
these, about 200 are singularity panels on the canard and wing surfaces. For 
such a panel scheme the required CPU time Is about 4.3 min on the CDC 7600. 
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DISCUSSION 


This study confirmed the validity of the approach adopted for wake-sheet 
modeling in that a good match with experimental data was obtained for the 
wing-alone case. Where a comparison was feasible POTFAN tended to give 
better overall performance than the pilot version of the APC. For the canard- 
wing case too, the results are encouraging. A more complete investigation 
should include body effects. 

It is postulated that improved results might be obtained by iterating on 
the leading-edge wake rollup (refs. 12 and 18) using the wake shapes generated 
by the Mangler-Smith theory only as an initial estimate. This would force 
the wake vortices to align more closely with the flow-field streamlines with 
successive iterations. Also to obtain better results for the pressure distri- 
butions, the Kutta condition at the trailing edge would have to be more 
accurate].y satisfied (refs. 19 and 20). This would necessitate rollup of the 
vortices from the rear edges as well as from the side edges. 


CONCLUSIONS 


The present investigation has demonstrated that the leading-edge wake 
estimate yields reasonable overall results for a modest investment of 
computer time. Accordingly, the method should be applicable to optimization 
schemes (ref. 21) to design canard-wing combinations. 
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(b) Doublet panel modeling. 
Figure 1.— Concluded. 
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(DISTANCES: cm) 



Figure 2.~ Placform view of closed coupled canard-wing configuration of refer- 
ence 13; ■ 4.12, A(. ■ 51.7“; iRy, ■ 2.5; ■ 60“ (distances: cm). 
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LIFT COEFFICIENT, Cl 



Figure 5.— Variation of lift with angle of attack for delta wing; yR ■ 1.15 

A ■ 74*; nominal panel density. 
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a POTFAN EXPERIMENT (REF. 15) 
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Figure 6.— Longitudinal distribution of normal force; ML ■ 1.15; A ■ 74*; 

nominal panel density. 
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LOCAL NORMAL FORCE COEFFICIENT 


POTFAN 

NO. OF CHORDWISE PANELS - 

16 

24 

32 



Figure 7.— Variation of longitudinal distribution of normal force with chord- 
wise panels; ffH ■ 1.15; A ■ 74*; number of spunwise panels ■ 8; a ■ 20®. 
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AXIAL VELOCITY AS FRACTION OF FREE-STREAM VELOCITY 



% LOCAL SEMISPAN 

(a) Axial component. 

Figure 9.— Spanwlse variation of velocity components above a wing at 
x/Cf ■ 0.8; Al ■ 1.07; A ■ 75*; nominal panel density a ■ 29". 
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SPANWISE VELOCITY AS FRACTION OF FREE-STREAM VELOCITY 



% LOCAL SEMISPAN 


(b) Spanwlse copponent. 
Figure 9.— Continued. 
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VERTICAL VELOCITY AS FRACTION OF FREE-STREAM VELOCITY 



(c) Vertical component. 
Figure 9.— Concluded. 
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TOTAL LIFT COEFFICIENT 



Figure 12.— Drag polar for canard-wing configuration; = A. 12; Ac = 51.7° 

= 2.5; A„ = 60°; M = 0.3. 
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□ EXPERIMENT (REF. 13) 
— POTFAN (WING + CANARD) 
APC (WING -t- CANARD) 



Figure 13.— Pltchlng-momeiit variation with angle of att/i'.k for canard-wing 
configuration; /Re - 4.12; Ac - 51.7"; - 2.5; A„ ^ 60"; M - 0.3. 
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PRESSURE COEFFICIENT ACp 



Figure 14.— Chordvise variation of pressure inboard of canard for canard-wing 
configurations with vortex rollup; - 4.12, Ag - 51.7“; ■ 2.5, 

A„ - 60“; M - 0.3; a - 28“. 
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rj - 0.8 (WITH RESPECT TO WING) 



Figure 15.— Chordwlse variation of pressure outboard of canard fur canard-wing 
configurations with vortex rollup; Ale ■ 4.12; ■ 51.7"; Aly ■ 2.5, 

Ay, - 60"; M - 0.3; a - 28". 
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